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ABSTRACT 

I Zwicky 1 (I Zw 1) is the prototype narrow line quasar. Its narrow line prohles 
minimize blending effects and thus allow a significantly more accurate study of individ¬ 
ual quasar emission lines. We report here the results of our study of the 1150-3250 A 
emission of I Zw 1 using a high S/N (50-120) spectrum obtained with the Faint Object 
Spectrograph (FOS) of the Hubble Space Telescope (HST). The UV spectrum of I Zw 1 
is very rich in emission lines, with essentially no wavelength range free of emission line 
features. The following main new results are obtained: 

1. The Mg II A2798 doublet is partially resolved, and the measured doublet ratio 
is 1.2/1. This ratio verifies theoretical predictions that the Mg II A2798 line is 
thermalized in the Broad Line Region (BLR). 

2. The A1 III A1857 doublet is clearly resolved, with a thermalized doublet ratio of 
1.25/1. The line optical depth provides an upper limit to the clouds distance from 
the ionizing source which is consistent with the “standard” BLR radius. 

3. A weak associated UV absorption system is detected in N V (possibly also in C IV 
and Lya as well), suggesting an outflow with a velocity of ~ 1870 km s“^ and 
velocity dispersion ^ 300 km s“^. 

4. Lines from ions of increasing ionization level show increasing excess blue wing 
flux, and an increasing line peak velocity shift, reaching a maximum blueshift of 
~ 2000 km s“^ for He II A1640. This may indicate an out-flowing component in 
the BLR, where the ionization level increases with velocity, and which is visible 
only in the approaching direction. The highest velocity part of this outflow may 
produce the associated UV absorption system. 

5. The small C III] A1909 equivalent width (EW), and the small C III] AlOOO/Lya 
and C III] A1909/Si III] A1892 flux ratios indicate a typical BLR density of Ue ~ 
10^^ cm“^, i.e. about an order of magnitude larger than implied by C III] A1909 
in most quasars. A BLR component with a density Ug 10^^ cm“^ is implied by 
the EW and doublet ratio of the A1 III A1857 doublet. 
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6. Prominent Fe II UV 191 emission is seen, together with weaker line emission at 
1294 A and 1871 A. These three features have been proposed by Johansson & 
Jordan (1984) as evidence for signihcant Lya pumping of the 8-10 eV levels of 
Fe II. 

7. Signihcant Fe III emission is present. The Fe III UV 34 and UV 48 multiplets are 
clearly resolved, and Fe III UV 1, UV 47, UV 50, and UV 68 may also be present. 
The implications of signihcant Fe III emission for the conditions in the BLR needs 
to be explored. 

8. Very weak [C III] A1907 and [Si III] A1883 emission may be present, suggesting a 
Narrow Line Region (NLR) component with Ug ~ 5 x 10^ cm“^. 

The rich UV emission spectrum of I Zw 1 should serve as a useful “template” for 
the identihcation of weak features in other Active Galactic Nuclei (AGNs), and as a 
useful benchmark for photoionization models, in particular for models of the complex 
Fe II emission spectrum. 


Subject headings: quasars:emission lines-quasars:absorption lines- quasars:individual 
(I Zwicky 1)- ultraviolet:galaxies 
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1. Introduction 

The prominent broad emission lines of quasars 
are typically blended with weaker emission fea¬ 
tures, such as emission from various Fe II multi- 
plets. This blending prevents a reliable study of 
the prominent line profiles, the measurement of 
weaker lines, and the identification of very weak 
lines. Even in high S/N spectra (e.g. Laor et al. 
1994a) the blending of features prevents a reliable 
study of the individual line profiles of even the 
most prominent lines. Blending effects are min¬ 
imized in the rare class of narrow line quasars 
(i.e. with FWHM;^ 1000 km s“^). Therefore the 
study of such objects can lead to a significantly 
more accurate study of the properties of individ¬ 
ual emission lines. 

I Zwicky 1 (I Zw 1) is the prototype narrow 
line Seyfert 1 galaxy (NLSl). Its optical spec¬ 
trum reveals narrow emission lines and strong 
Fe II emission (Sargent 1968; Phillips 1976, 1978; 
Oke & Lauer 1979). The NLSl class of AGNs 
also tend to have weak [O III] AA4959, 5007 emis¬ 
sion, asymmetric H/3 profiles (Boroson & Green 
1992, hereafter BG92), steep soft X-ray spec¬ 
tra (Laor et al. 1994b, 1997a; Boiler, Brandt 
& Fink 1996), and possibly strong IR emission 
(e.g. Lipari 1994). Although I Zw 1 is classi¬ 
fied as a Seyfert 1 galaxy, its absolute luminos¬ 
ity, My = —23.8 (for Hq = 50 km s“^ Mpc“^, 
( 7 o = 0), actually qualifies it as a low luminosity 
quasar (Veron-Cety & Veron 1991), we therefore 
refer to it as a quasar throughout the paper. 

In this paper we describe a detailed study of 
the UV emission line spectrum of I Zw 1 obtained 
with the Faint Object Spectrograph (FOS) of the 
Hubble Space Telescope (HST). An earlier UV 
spectrum of I Zw 1 obtained with the Interna¬ 
tional UV Explorer is described by Wills (1983), 
and HST UV spectropolarimetry of I Zw 1 is de¬ 
scribed by Smith et al. (1997). Our study has 
two main goals. The first goal is to provide a de¬ 
tailed “identification template” of the UV emis¬ 
sion features present in the spectra of quasars. 


The second goal of this study is to provide a 
high S/N ratio Fe II UV emission template. The 
strong and narrow Fe II emission lines in I Zw 1 
make it an ideal object for the construction of 
such a template, as was clearly demonstrated in 
the optical regime by BG92. Such a template is 
also important for studies of UV Fe II emission in 
AGNs (e.g. Netzer k, Wills 1983; Wills, Netzer, 
k Wills 1985; Collin-Souffrin k Dumont 1986). 
These studies are important for understanding 
the physics of the broad-line region (BLR) since 
Fe II is a major, and sometimes the dominant 
coolant of the BLR (e.g. Phillips 1977, 1978; 
Miley k Miller 1979; Bergeron k Kunth 1984; 
Joly 1981; Lipari, Terlevich, k Macchetto 1993). 
The Fe II template is also important for stud¬ 
ies of lines strongly blended with Fe II, such as 
Mg II A2798. 

In this paper we concentrate on the first goal 
of this project. In §2 we describe the observa¬ 
tions and present the “identification template”. 
In §3 we discuss various diagnostics provided by 
some of the UV lines, and in §4 we describe the 
physical properties of an associated UV absorp¬ 
tion system, and discuss its possible relationship 
to observed UV emission and X-ray absorption. 
A brief discussion of the similarities of I Zw 1 to 
broad absorption line quasars (BALs) (§5) and a 
summary of our main conclusions (§6) complete 
the paper. 

2. Observations and Analysis 

We obtained high S/N spectra (50—120 per 
resolution element) of I Zw 1 with the three high 
resolution (R=1300) gratings of the FOS from 
1150 A to 3280 A. This S/N is comparable to the 
highest S/N yet obtained with HST for bright 
AGNs. All the observations were made using 
the 0”.86 aperture. The dates of the observa¬ 
tions are; G130H, 13 February 1994 (UT Date); 
G190H and G270H, 14 September 1994. The 
HST spectra were reduced as described in Schnei¬ 
der et al. (1993). The wavelength zero points of 
the three HST spectra were set using Galactic 
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interstellar medium absorption lines (Savage et 
al. 1993). The size of the necessary shift and the 
lines used for each of the three spectra (G130H, 
G190H, and G270H) were respectively 0.36 A 
(Si II A1526.71), 1.08 A (Al II A1670.79), and 
0.99 A (Fe II lines). We also obtained a compli¬ 
mentary high resolution ground based spectrum 
which covers the range from 3183 A to 4074 A. 
The spectrum was obtained using the KPNO 
2.1m telescope and Goldcam spectrograph on 22 
September 1995 UT. The grating used was a 
837 lines/mm grating, used in 2nd order, pro¬ 
viding a dispersion of approximately 0.47 A per 
pixel and a spectral resolution of 1.2 A. The 
Ca II AA3933.66, 3968.47 doublet absorption by 
the ISM was used to secure a zero point wave¬ 
length which matches the UV zero points, which 
were also set using lines of low ionization state 
ions. The optical spectrum is used in order to 
get complete coverage of the 2000-4000 A region, 
required for generating the Fe II template. 

Figure 1 presents the overall optical-UV spec¬ 
trum of I Zw 1 from llOOA to 6000A. The spec¬ 
trum was corrected for the effects of Galactic 
extinction of E(B—V)=0.1 using the standard 
Seaton (1979) reddening curve. The extinction 
is deduced from the neutral hydrogen column 
density of Nh i = 5.07 x 10^^ cm“^ taken from 
Elvis, Lockman, & Wilkes 1989), and the relation 
E(B— V)=Nh/4.93 X 10^^ found by Diplas & Sav¬ 
age (1994). The rest wavelength was calculated 
assuming z = 0.0608, the redshifts of Ha and 
H/I (e.g. Phillips 1976). A composite spectrum 
of ~ 700 quasars observed from the ground by 
Prancis et al. (1991, hereafter the Erancis et al. 
composite) is also shown for comparison. This 
composite represents well the “average” quasar’s 
continuum, and matches well similar composites 
constructed for lower z quasars with HST (see 
Pig.2 in Laor et al. 1995). 

The optical spectrum of I Zw 1 displayed in 
Pig.l was kindly provided by B. Wills (1996, pri¬ 
vate communication). It matches the HST spec¬ 
trum in the overlap region, 3170-3270A, to bet¬ 


ter than 5%, and is used here to demonstrate 
the overall spectral shape. A flux discontinuity 
may be present just blueward of G IV (logzz = 
15.29) since the region at logiy > 15.29 (ob¬ 
tained with the G190H) was observed 7 months 
after the logiy < 15.29 region (obtained with the 
G130H). The overall optical-UV power-law slope 
is a ~ — 1, significantly redder than the slope of 
the typical quasar (a —0.3, as displayed by 
the Prancis et al. composite or other optically 
selected samples, e.g. Neugebauer et al. 1987). 
This may indicate that I Zw 1 is significantly 
reddened, as further discussed in §3.3. Note that 
many weak features in the Prancis et al. compos¬ 
ite are seen in I Zw 1 as clear and well resolved 
lines. 

Figures 2a-2c present an expanded view of the 
UV spectrum shown in Fig.l. The rest frame 
positions of various lines are marked above the 
spectrum. Lines marked below the spectrum 
designate various absorption lines originating in 
the Galactic interstellar medium (ISM). All ISM 
lines detected here are typically seen in FOS 
spectra of quasars (Savage et al. 1993). We make 
no attempt here to identify the many Fe H emis¬ 
sion blends present. The various Fe H multiplets 
designated in Fig.2 are taken from Moore (1952), 
and serve only as an illustration of possible Fe H 
features. The number of possible Fe H multiplets 
is extremely large, and a reliable identification 
requires detailed theoretical calculations, which 
are beyond the scope of this paper. 

Figure 2d presents the optical spectrum. The 
spectrum was obtained in a single exposure and 
is affected by cosmic ray hits, which were deleted 
by taking out 5 small patches (width 1-2 A) from 
the spectrum. The absolute flux calibration is 
not secure. The wavelength zero point was set us¬ 
ing Ga H AA3934.78, 3969.59 Galactic interstellar 
medium absorption lines. Since the HST wave¬ 
length zero point is also based on low ionization 
species (Si H, Al H, Fe H), it probably matches 
well the optical spectrum zero point. The HST 
FOS G270H spectrum which overlaps with the 
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optical spectrum in the 3183 A to 3280 A re¬ 
gion is also shown in Fig.2d, shifted down by a 
factor of 0.37 in log flux. Apart from the flux dis¬ 
crepancy, the spectral features match very well. 
All the optical Fe II multiplets from the first 20 
optical Fe II multiplets listed by Moore (1959) 
are marked for illustration only (it does appear 
though that some of the observed lines agree with 
expected positions). Note the absence of the 
[O II] A3727 and [Ne V] AA3346, 3426 forbidden 
lines . 

Table 1 lists for each emission feature the mea¬ 
sured rest frame wavelength (defined as the wave¬ 
length of the peak flux in the Gaussian smoothed 
spectrum), the expected wavelength, the emit¬ 
ting ion, the observed flux, rest frame equiva¬ 
lent width (EW), full width at half maximum 
(FWHM), and the velocity shift of the peak of 
the emission line. These parameters were mea¬ 
sured by fitting a set of Gaussians to each emis¬ 
sion feature, as described in Laor et al. (1994a, 
1995). In strongly blended features, or in very 
weak features, some of the emission parameters 
could not be measured reliably, and are therefore 
not listed. We preferred not to assign errors to 
the measured parameters since such errors may 
not be very meaningful given the very compli¬ 
cated overall spectrum and the subjectivity of 
both the placement of the continuum and sepa¬ 
ration of blended features. 

Figure 3 presents a comparison of an Ha pro¬ 
file to the observed profile of most of the strong 
emission features. Direct comparison of line pro¬ 
files is not meaningful in most cases since most 
lines are strongly blended (despite their narrow 
profile in I Zw I). The comparison is therefore 
made using a template profile, which allows us to 
test if the observed profile can be reconstructed 
with a sum of template components represent¬ 
ing the various lines which may contribute to the 
blend. For the template we use the high S/N 
spectrum of Ha obtained for I Zw 1 as part of 
the BG92 study. The Ha line is used since it 
has the highest S/N, and is least blended. The 


wavelength used to set the velocity scale is given 
in each sub-panel. Most features are blends, and 
in each subpanel we list the ratio of fluxes in 
the blend components assumed in the fit. When 
searching for the best fit we allowed the con¬ 
tinuum level to be below the apparent contin¬ 
uum level, since the observed continuum may be 
contaminated by weak unresolved emission (see 
in particular the Mg H doublet fit, discussed in 
§3.1.5). Note that the Ha profile is slightly asym¬ 
metric with some excess flux on the blue wing 
(e.g. see the C IV 1/1 doublet fit). The tem¬ 
plate fits were not used for measuring any of the 
emission line parameters. 

3. Emission Lines 

Below we describe some of the BLR diagnos¬ 
tics provided by ratios of lines from different lev¬ 
els of a given ion, or within certain multiplets of 
a given ion. 

3.1. Line Diagnostics 

3.1.1. He 

The Rueii =He H A4686/He H A1640 flux ra¬ 
tio can be used as a reddening indicator. Below 
we explore its usefulness in I Zw 1. There is a 
very weak, if any, feature at the expected He H 
rest wavelength of 1640.4 A. The nearest strong 
feature is at 1629.5 A, which we suggest is highly 
blueshifted ( 2000 km s ^) He H emission. Al¬ 

ternatively, this feature may be just another Fe H 
feature, many of which are scattered throughout 
the spectrum. However, we And the He H iden¬ 
tification more likely since C IV and N V also 
display large blue shifts (~ 900 km s“^), and the 
He H velocity shift is typically twice that of G IV 
(e.g. Laor et al. 1995, Table 7). In addition the 
He H/Lya flux ratio of ~ 0.03 measured using 
the 1629.5 A feature is typical of quasars (e.g. 
Laor et al. 1995, Table 6). 

Photoionization models indicate values of about 
7?Heii — 8 — 11 (Netzer , Elitzur, & Ferland 1985). 
BG92 did not detect the He H A4686 line. A 
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plausible upper limit on the EW is 0.5 A using 
the BG92 high S/N optical spectrum, which im¬ 
plies a line flux< 5 x 10^® erg s“^ cm“^ A“^, 
where the continuum flux level is obtained from 
a lower resolution optical spectrum presented in 
Fig.l. We thus get i?Heii > 18. However, since 
the He H A1640 line is very broad and asymmet¬ 
ric, the He H A4686 upper limit may have been 
somewhat underestimated, while some of the flux 
attributed to He H A1640 may be due to blended 
Fe H features. Therefore, the i?Heii does not pro¬ 
vide a useful contraint on the reddening. 

3.1.2. C 

Metal lines are produced mostly by collisional 
excitation of ions in their ground level. The 
C HI* A1176 line discussed below is unique in 
originating from radiative decay between two rel¬ 
atively high lying levels. The transition is from 
the 2p^ level, 17 eV above the ground level, 
down to the metastable 2s2p^P level, 6.5 eV 
above the ground level (which decays to ground 
level via the C HI] A1909 line). The typical tem¬ 
perature in the C HI region in photoionized gas 
is < 2 X 10^ K, and the small Boltzman factor, 
g-x/fc'T < 5 X 10“® (x = 17 eV), suggests that 
the collisional excitation rate would be low and 
therefore that the line flux should be small. 

The C HI* A1176 line is clearly detected in 
I Zw 1. The presence of this line was already sug¬ 
gested by Laor et al. (1995, Fig.4) in high S/N 
FOS spectra of PG 1116-1-215, PG 1216-1-069, 
and PG 1538-1-477, but these detections were 
only marginal due to the effects of absorption 
by the Lya forest, and due to strong blending 
with the Lya blue wing. These two effects are 
negligible in I Zw 1, and the line can therefore 
be clearly detected. The C HI* A1176/Lya flux 
ratio found in those three objects is 0.02 — 0.05 
versus a comparable ratio of 0.014 found here. 
This line was also detected recently in a compos¬ 
ite HST spectrum of quasars at a similar level 
(Hamann et al. 1997). 

Laor et al. (1997b) briefly discussed vari¬ 


ous excitation mechanisms for C HI* A1176, and 
showed that the observed high C HI* A1176/ 
G HI] A1909 flux ratio is ~ 50 times larger than 
expected if collisional excitation at typical BLR 
conditions is the only significant process. How¬ 
ever, as discussed below, the BLR density in the 
G HI region of I Zw 1 is probably about an order 
of magnitude larger than in a typical BLR, which 
strongly suppresses C HI] A1909 emission. In ad¬ 
dition, the gas temperature in the BLR tends 
to increase with increasing density (due to ther- 
malization of some of the main coolants), which 
increases the C HI* A1176 emission, as shown 
by Korista et al. (1997, their Figure 3c). The 
suppression of C HI] A1909 and enhancment of 
G HI* A1176 at n-e 10^^ cm“^ can explain 
the high C HI* A1176/G HI] A1909 flux ratio 
in I Zw 1. Note that dielectronic recombination 
(Storey 1991) also contributes significantly to the 
G HI* A1176 flux. Laor et al. (1997b) suggested 
that dielectronic recombination is not significant 
in I Zw 1 based on the C HI A2297 line flux, 
but this suggestion was based on a misquote of 
Storey (1991), and the actual values do not rule 
out a significant contribution from dielectronic 
recombination. Resonance scattering of contin¬ 
uum photons may also contribute significantly 
to the C HI* A1I76 emission if the velocity dis¬ 
persion within the BLR clouds is large enough 
(Ferguson, Ferland & Pradhan 1995; Laor et al. 
1997b). 

The [C HI] A1907 line may be present in 
I Zw 1. This line is produced by a forbidden tran¬ 
sition from the ^P 2 term in the 2s2p conhguration 
of C HI down to the ground level, with a critical 
density of Uc = 5 x 10^ cm“^ (Osterbrock 1989, 
Table 3.11). The ^P 2 term lies 0.007 eV above 
the ^Pi term which produces the well known 
G HI] A1909 line, where ric ^ 5x 10® cm“^. Thus, 
the flux ratio Rq hi =[C HI] A1907/G HI] A1909 
provides a temperature independent density di¬ 
agnostic (e.g. Osterbrock 1989, Fig. 5.5). 

The observed C HI] A1909 line profile appears 
to be double peaked with the short wavelength 
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peak at 1906.4 A matching well the expected 
1906.68 A position of the [C III] line. This iden¬ 
tification is not secure since excess emission is 
seen on the blue side of some of the other emis¬ 
sion lines (see §3.2), and the significance of the 
feature is only marginal. But, an analogous [Si 
III] A1883 line may also be present (see §3.1.7), 
which implies similar conditions. If true, then the 
[C III] line is much narrower and weaker than the 
C III] line. The observed Rq hi ^0.1 then im¬ 
plies that some of the C III emission originates in 
gas with an electron density ng ^ 5 x 10^ cm“^. 

The only clear detection of the [C III] A1907 
line in AGNs is apparent in the high-dispersion 
lUE spectrum of NGC 4151 displayed by Bro- 
mage et al. (1985, Fig.l). 

The G III] A1909 line is very weak in I Zw 1. 
Both its EW of ~ 2.5 A, and the G III] A1909/ 
Lya=0.023 flux ratio are a factor of 5-10 smaller 
than typically observed (e.g. Francis et al. 1991; 
Laor et al. 1995; Wills et al. 1995). Since 
C III] A1909 is thermalized at rig 5 x 10® cm“^, 
the suppression of C III] A1909 may imply that 
the typical BLR density in I Zw 1 is about 
an order of magnitude larger than indicated by 
C III] A1909 in typical AGNs, i.e. rig ~ 10^^ cm“^ 
instead of ng ~ 10^® cm“^. Additional evi¬ 
dence for this interpretation is provided by the 
Si III] A1892 discussed in §3.1.7. 

3.1.3. N 

The N V AA1239,1243 doublet components 
are separated by the equivalent of ~ 1000 km s“^, 
and the narrow line profiles in I Zw 1 suggest that 
this doublet should be resolvable. However, the 
N V profile shows a large excess flux at the blue 
wing, and the Ha template fit is inadequate (see 
Fig.3). It does appear, however, that the doublet 
flux ratio is significantly lower than the optically 
thin 2/1 ratio. The physical implication of the 
observed doublet ratio is discussed below for the 
Mg H and Al HI doublets which are rather well 
resolved (§§3.1.5, 3.1.6). 


3.1.4. O 

The O I A1304 line is thought to be produced 
through Bowen fluorescence. In this process a 
Ly/3 (1025.72 A) photon is absorbed by O I in 
the ground 2p^P level, and gets excited to the 
3d^D level, which is within the thermal core of 
Ly/3 (1025.77 A). The decay back to the ground 
level is either direct, or through two intermedi¬ 
ate levels. The indirect decay route results in 
a A11287 photon {3d^D to 3p^P), a A8446 pho¬ 
ton {3p^P to 3s^S), and a A1304 photon (3s^5 
to 2p^P). Thus, one expects equal number of 
photons in these three lines. Rudy, Rossano & 
Puetter (1989) very nicely confirmed this predic¬ 
tion for the A11287 and A8446 lines in I Zw 1. 
Below we show that the same number of photons 
is not seen in the A1304 line, and briefly discuss 
possible implications. 

The O I A1304 line is blended with the Si H AA 
1304, 1309 doublet. This blending is very strong 
in most quasars preventing a reliable estimate 
of the relative contributions of O I and Si H to 
the blend even in studies of high quality spectra 
dedicated to Si H (e.g. Appendix A in Baldwin 
et al. 1996). In I Zw 1 the Si H A1309 com¬ 
ponent of the Si H doublet is clearly resolved. 
The Si H A1304 component is hopelessly blended 
with the O I A1304 line, which is actually a triplet 
at 1302.17, 1304.86, and 1306.0 A. However, the 
Si H 1309/1304 flux ratio should either be 2/1 
(optically thin case), or 1/1 (optically thick case), 
and thus the clear detection of the 1309 compo¬ 
nent allows an estimate of the 1304 component 
flux within a factor of two. The template fit to 
the O I-|-Si H blend (Fig.3), which assumes an 
optically thick Si H doublet, suggests that 0.5 of 
the blend flux is due to O I. If the Si H doublet 
is optically thin (less likely given the results in 
§§3.1.5, 3.1.6), then the Si H A1304 component 
is only half as strong, and then 0.56 of the blend 
flux is due to O I. Using the O I A8446 flux from 
Rudy et al. then gives Rq i = A8446/A1304 pho¬ 
ton flux ratio is 5.5 (4.8), for the optically thick 
(thin) case. 
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The shortest wavelength flux in the spectrum 
from which the Rudy et al. O I A8446 flux is 
taken (at 7500A) agrees within ~ 20% of the 
flux in the optical spectrum shown in Fig.l, thus 
variability between the times the A8446 and the 
A1304 lines were measured is most likely not a 
significant source of error. 

Kwan & Krolik (1981) have found that Rq i 
can increase from unity to 1.3 due to Balmer 
continuum absorption of A1304 photons, and the 
production of A8446 photons by collisional ex¬ 
citation of O I ions in the 3s level. A third 
mechanism which can increase Rq i, noted by 
Grandi (1983), is the significant probability for 
a 2p^^D - 3s^S A1641 semi-forbidden transition 
for rAi304 > 10^, or equivalently TLya > 10®, 
which would further increase Rq i- Kwan (1984) 
included the O I] A1641 transition and indeed 
found Rq I ~ 3.2, where most of the increase 
from unity is due to the “re-routing” of A1304 
photons to A1641. 

This mechanism can be tested in I Zw 1 since 
the He II A1640 line is strongly blueshifted (to ~ 
1630 A), while the O I lines are not. This allows 
us to look for the O I] A1641 line and test if it can 
indeed explain the very high Rq i observed. No 
strong line is detected at 1641 A, and we estimate 
that such a line cannot add more than ~ 30% 
to the observed A1304 line. The absence of a 
significant flux in the O I] A1641 line therefore 
implies that most of the Ly/3 absorption by O I 
in I Zw 1 takes place at a depth where TLya < 10®. 

Rudy et al. gave observational evidence that 
collisional excitations cannot significantly increase 
Rq I in I Zw 1. Thus, of the three mechanisms 
mentioned above, the only one left is Balmer con¬ 
tinuum absorption of A1304 photons. It remains 
to be explored whether this mechanism can sup¬ 
press the A1304 line by the required factor of ~ 5 
when TLya < 10®. 

Another option is that there are few O I A1304 
photons because the observed emission is strongly 
reddened, as discussed below in §3.3. 


Another observed emission line from O is the 
O IV] A1402 line. This line is a multiple! of five 
components and it is strongly blended with the 
Si IV A1397 doublet producing a single broad 
A1400 feature in most AGNs. The relative con¬ 
tribution of O IV and Si IV has been debated in 
the past (Baldwin & Netzer 1978; Wills & Netzer 
1979; Tytler & Fan 1992; Wills et al. 1993; Laor 
et al. 1994a). The narrow lines of I Zw 1 ap¬ 
pear to provide a good opportunity to determine 
accurately the relative contribution of O IV and 
Si IV to the A1400 feature. 

The peak of the A1400 blend is clearly split in 
I Zw 1 into two peaks at 1392.1 A and 1400.6 A. 
These correspond well to the expected positions 
of the Si IV doublet at 1393.76 A and 1402.77 A, 
blueshifted by ~ 2 A. However, some O IV] must 
also be present based on the following argument. 
The 1392.1 A feature must be due to Si IV only 
(the shortest wavelength O IV] component is at 
1397.2 A). The second Si IV doublet component 
is at 1400.6 A and its flux should either be equal 
to (optically thick) or half of (optically thin) the 
1397.2 A line flux. The observed 1400.6 A com¬ 
ponent is stronger than the 1392.1 A component, 
and must therefore include some O IV] emis¬ 
sion. A specific template deblending is shown 
in Fig.3, using Ha as a template, which implies a 
Si IV/0 IV] flux ratio of 1/2. However, this ra¬ 
tio is not well constrained since the intrinsic line 
profile is probably not symmetric, as suggested 
by the large observed blue excess flux, and the 
fact that all other high ionization lines display a 
strong blue excess flux (see §3.2). Thus, we can 
only verify that both Si IV and O IV] are present, 
and that they have roughly comparable flux. 

3.1.5. Mg 

The Mg H ion is produced by photons above 
7.646 eV and destroyed by photons above 15.035 eV, 
it therefore survives mostly in the partially ion¬ 
ized H region, where photons above 13.6 eV are 
blocked. The Mg H A2798 doublet is therefore a 
diagnostic of the conditions in the extended par- 



tially ionized layers inside the BLR clouds. The 
flux ratio of the Mg II 2796.35 A/2803.53 A com¬ 
ponents indicates the optical depth in the line. 
Below we describe the observed flux ratio and 
discuss its implication. 

The Mg II A2798 profile is clearly double 
peaked. Its core is well fit with a 1.2/1 ratio for 
the 2796.35 A/2803.53 A doublet lines using the 
Ha profile as a template (Fig. 3). A good fit re¬ 
quires that the actual continuum level be placed 
20% below the observed flux level at the appar¬ 
ent line base. This continuum level matches very 
well the continuum level inferred at 2800 A from 
a linear interpolation of the continuum at 2650 A 
and 3020 A. 

What does the 1.2/1 doublet ratio implies for 
the conditions in the Mg II emitting gas? The 
2803.53 A, 2796.35 A lines correspond to tran¬ 
sitions from the ‘ip^Pi /2 and 3p^P3/2 levels to 
the ground 35^51/2 level. The collision strength 
at T = 10^ K is H(2p,2 5) = 16.9 (Pradhan 
& Peng 1995), yielding collisional de-excitation 
rates of q{3p‘^P3/27^s'^Si/2) = 3.65 x 10“^ and 
q{Sp‘^Pi/ 27 ^s‘^Si/ 2 ) = 7.3 X 10“’^. The radiative 
de-excitation rates for the ^p‘^P 3 i 2 ,il 2 levels are 
A = 2.6 X 10® s“^ (Morton 1991). Collisional 
de-excitation dominates when 

Tone > A/q, 

where tq is the line center optical depth for reso¬ 
nance scattering, i.e when ToUe > 7.1 x 10^^ cm“® 
and 3.6 x 10^^ cm“® for the 3p2pg^2 and ^p‘^Pi /2 
levels. For Tone > 7.1 x 10^^ cm“® the Mg II 
2796.35 A/2803.53 A doublet ratio is thermal- 
ized, i.e. it changes from 2/1 to 1/1 (see Hamann 
et al. 1995 for results on other analogous UV 
doublets). Thus, the observed 1.2/1 ratio implies 
that most of the Mg II emission is produced in a 
region where Tong ^ 10^® cm“®. 

What does the above result implies for the 
conditions in the partially ionized BLR gas? As 
shown below, the value of Tone provides a di¬ 
rect measure of the distance of the BLR clouds 
from the central ionizing continuum source. The 


line center optical depth is tq = (Ti,pSMg ii where 
= 1.47 X 10“^2 (3)332 -g total Mg II dou¬ 
blet line center absorption cross section (for T = 
10^ K), and ^Mg II is the column density of Mg II 
ions. The partially ionized H region, where Mg II 
survives, has a column density about 10-30 times 
larger than the column density of the highly ion¬ 
ized H region, given by Sh ii = lO^®!/ cm“2^ 
where U is the ionization parameter (= n^jne, 
the photon/electron density). Balmer contin¬ 
uum absorption may limit the depth of the Mg II 
doublet emitting layer, but this absorption will 
not affect our conclusion. Using the Mg/H solar 
abundance ratio of 3.8x 10“®, and assuming most 
Mg is Mg II, we get ^Mg n ~ 4 x 10^^U 
or To = 6 X 10^17, assuming the gas is ionization, 
rather than matter bounded. The Mg II doublet 
is therefore thermalized when 6 x lO^I/ne ^ 10^®, 
or ^ lO'^ cm“®. 

Reverberation mapping (e.g. Netzer 1990; Pe¬ 
terson 1993; Maoz 1995; Kaspi et al. 1996), 
and theoretical considerations (Netzer & Laor 
1993) suggest Rblr = O-lL^g pc, which implies 
~ 2 X10® cm“®in the BLR. Thus, the observed 
1.2/1 Mg II doublet ratio provides direct confir¬ 
mation that the BLR is indeed optically thick 
in the Mg II A2798 line, as predicted by ‘stan¬ 
dard’ photoionization models for the BLR. More 
refined calculations using actual photoionization 
models can be used to explore the exact doublet 
ratio predicted by various BLR models and its 
agreement with the observed value. 

Note that one generally expects a doublet ra¬ 
tio from optically thick gas to be somewhere be¬ 
tween 2/1 and 1/1, rather than 1/1, since (as 
is always the case in atmospheres), most of the 
contribution to the line flux arises from the re¬ 
gion where the emission just becomes optically 
thick, i.e. from the region where the doublet ra¬ 
tio changes from 2/1 to 1/1. 

3.1.6. Al 

The components of the Al III A1857 doublet 
are widely spaced (~ 1300 km s“^) and are there- 
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fore clearly resolved in I Zw 1. A1 III is created 
by photons above 18.829 eV, and destroyed by 
photons above 28.448 eV, thus unlike Mg II, it 
exists in the highly ionized H region only, and the 
observed A1 III doublet ratio provides constraints 
on the BLR H II region, as discussed below. 

The template ht suggests a 1.25/1 ratio for 
the 1854.71 A/1862.79 A doublet lines (Fig.3). 
This ratio is somewhat uncertain since there are 
significant deviations from the Ha template fit, 
possibly due to blending with Fe II multiplets 
(e.g. Wills, Netzer &: Wills 1980), or an intrinsic 
blue excess flux, as displayed by other high ion¬ 
ization lines (§3.2). Baldwin et al. (1996) found 
a similar 1/1 flux ratio in three objects where the 
Al III doublet was clearly detected. 

Below we calculate the expected Al III A1857 
doublet ratio, as calculated above for the Mg II 
A2798 doublet. 

The Al III and Mg II ions are both Na I 
like. Thus, the 1854.71 A and 1862.79 A Al III 
lines are also produced by transitions from the 
3p^Pi/2 3p^P3/2 levels to the ground 3s ^5^/2 

level. The collision strength at T = 10^ K is 
S) = 16.0 (Dufton & Kingston 1987), the 
radiative de-excitation rate is A = 5.4 x 10® s“^ 
(Morton 1991), the solar Al/H abundance ra¬ 
tio is 3.6 X 10“®, and a,yQ = 9.4 x 10“^® cm^. 
Al III exists in the highly ionized H region only, 
where Sh ii = 10^® t/ cm“^. Using these param¬ 
eters, and assuming Al is mostly Al III, we find 
that the Al III A1857 doublet is thermalized for 
n..y ^ 4 X 10® cm 

The observed 1.25/1 Al III doublet ratio sug¬ 
gests that Al III is formed in clouds at, or inside 
the ‘standard’ BLR radius of 0.1L4g pc, where 
^ 2 X 10® cm“®. The fact that the Al III 
doublet is expected to be only marginally ther¬ 
malized at the ‘standard’ BLR radius suggests 
that this doublet ratio can provide an interesting 
probe of the size of the BLR. 

The relatively large Al III A1857 EW (4.4 A) 
requires a BLR component where Al III is the 


dominant ionization state of Al, which occurs for 
U ~ 10“^ — 10“® (e.g. Baldwin et al. 1996, 
Fig.7a). Since ^ 4 x 10® cm“®, the above 
range in U implies Ug ^ 4 x 10^^ — 4 x 10^^ cm“®, 
which is comparable, or denser, than the com¬ 
ponent responsible for the C III] A1909 emission 
(§3.1.2). 

We caution that resonance scattering of con¬ 
tinuum photons has not been included in the 
above discussion. This process could contribute 
an EW of ~ 2 A to each of the Al III A1857 dou¬ 
blet components, i.e. as much as observed, if the 
velocity dispersion in the BLR clouds is as large 
as the lines FWHM, and the clouds cover ^ 30% 
of the continuum source (e.g. Eerguson, Ferland 
& Pradhan 1995). Continuum resonance scatter¬ 
ing would also imply a much lower value for tq 
than estimated above (due to the large velocity 
gradients required in the BLR clouds). Calcu¬ 
lations with detailed photoionization codes are 
required in order to understand whether contin¬ 
uum resonance scattering is indeed viable. 

3.1.7. Si 

The Si II 1260.42 A, 1264.74 A, 1265.00 A, 
1304.37 A, 1309.28 A, 1808.01 A, 1816.93 A, and 
1817.45 A lines are clearly detected. The flux 
of the Si II A1194 and A2335 blends is well con¬ 
strained. This wealth of data provides a unique 
opportunity to explore the Si II line formation 
mechanisms (e.g. Baldwin et al. 1996, Appendix 
C). Such analysis is beyond the scope of this pa¬ 
per. 

The outer shell electronic configuration of the 
Si III ion is analogous to that of C III. The anal¬ 
ogous line to the [C III] A1907 line (§3.1.2) is the 
[Si III] A1883 line, which can also be used as a 
density diagnostic for the Si III gas. Weak [Si 
III] A1883 emission may be present at 1881.1 A. 
The identification of this marginal feature is sup¬ 
ported by the fact that its 250 km s“^ blue shift 
matches the Si III] A1892 shift, and by the possi¬ 
ble presence of the analogous [C III] A1907 line. 
As for the C III ion, the Rsi m = [Si III] A1883/ 
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Si III] A1892 flux ratio is a robust density indi¬ 
cator (Keenan, Feibelman & Berrington 1992). 
The observed i?si iii ^ 0.1 implies that some of 
the Si III resides in gas with Ue ~ 5 x 10® cm“®, 
i.e. the same conditions implied by Rc iii, as ex¬ 
pected since the spatial distribution of both ions 
in photoionized gas should largely overlap. 

The Si III] A1892/C III] A1909 flux ratio of 
~ 3.5 is significantly larger than the typical ra¬ 
tio of ~ 0.3 ±0.1 observed in quasars (Laor et 
al. 1995). This high ratio results from the fac¬ 
tor of 5-10 suppression in the C III] AI909 flux 
(§3.1.2), rather than significant enhancement of 
Si III] A1892. The most likely interpretation for 
this high ratio is a relatively dense BLR. The crit¬ 
ical density for C III] A1909, n-c ~ 5x 10® cm“®, is 
significantly smaller than for Si III] A1892, where 
Uc = 1.1 X 10^^ cm“^. Thus, if the BLR density is 
10^^ cm“®, rather than 10^® cm“^, C III] A1909 
would be suppressed by a factor of ~ 10, while 
Si III] A1892 would not be affected, consistent 
with our observations. 

3.1.8. Fe 

The positions of some Fe II UV blends are 
marked in Fig.l only to illustrate various possi¬ 
ble Fe II emission features. Numerous additional 
weak Fe II blends are most likely present at other 
wavelengths. Due to the extreme complexity of 
the Fe II ion we do not attempt to identify indi¬ 
vidual features. 

The likely presence of Fe II emission well be¬ 
low 2000 A, and possibly even at ~ 1110-1130 A, 
indicates that Fe II is excited by processes other 
than just collisions, as the electron temperature 
in the Fe II region is most likely far too low for 
significant collisional excitation of levels 10 eV 
above the ground state. One such process, sug¬ 
gested by Penston (1987), is resonance absorp¬ 
tion of Lya photons by Fe 11. Johansson & Jor¬ 
dan (1984) found Lya resonance absorption to be 
significant in various stellar systems, and identi¬ 
fied the Fe II A1294, Fe II UV 191 A1787, and 
the Fe II A187I multiplets as the signatures of 


such a process (see also Johansson &: Hamann 
1993). The Fe II UV 191 multiplet is very promi¬ 
nent in I Zw 1, and weak emission blends are 
clearly apparent at 1294 A, and 1870 A, sug¬ 
gesting that resonance absorption of Lya pho¬ 
tons may be a significant excitation mechanism 
for Fe II in I Zw 1 as well. Resonance scatter¬ 
ing of continuum photons, and Fe II-Fe II line 
flouresence could also be a significant process for 
populating high lying Fe II levels (e.g. Netzer & 
Wills 1983; Wills, Netzer & Wills 1985). Clearly, 
the Fe II rich spectrum of I Zw 1 should serve as 
a valuable tool for future studies of Fe II emission 
in AGNs. 

A number of Fe III multiplets are clearly 
present. The Fe III UV 34 multiplet, clearly 
identified here, was first discovered by Hartig 
& Baldwin (1986) in a broad absorption line 
quasar. An emission feature near 2070 A 
is commonly seen in quasars (e.g. Wills et al. 
1980), and was identified as a likely Fe II blend. 
Baldwin et al. (1996) suggested that this fea¬ 
ture is due to Fe III UV 48 emission. Here the 
three Fe III UV 48 components at 2062.2 A, 
2068.9 A, and 2079.65 A are clearly resolved, 
which verifies the identity of this feature. The 
Fe III UV 47 AA2419.3,2438.9 doublet is most 
likely present as well. Other Fe III blends, such 
as Fe III UV 50 and Fe III UV 68, and in particu¬ 
lar the resonance Fe III UV 1 at 1122-1130 A may 
also be present, but these identifications cannot 
be verified here since these blends are not clearly 
resolved. 

The relative flux ratio of the multiplet com¬ 
ponents can provide an important diagnostic for 
the optical depth and excitation mechanism of 
Fe III. Given the complexity of the Fe III ion, 
these diagnostics are beyond the scope of this 
paper. 

3.2. Line Profiles and Velocity Shifts 

Optical spectra of I Zw 1 revealed two velocity 
systems. The first system is at z = 0.0608 for the 
low ionization forbidden lines and Balmer lines, 
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including Ha, H/3, optical Fe II multiplets, [Ca II] 
AA7291.46, 7323.89, [S II] AA4068.60,4076.35, 
[N II] A5754.57, He I A5875.6, Na I AA5889.95, 
5895.92, and [O I] AA6300.3,6363.8. The sec¬ 
ond system is at z = 0.0587 (i.e. blueshifted by 
~ 630 km s“^) for the higher ionization forbid¬ 
den lines, [O HI] AA4958.82, 5006.85 and [Ne HI] 
AA3868.74,3967.47 (Phillips 1976, Oke & Lauer 
1979). 

We find the same trend in the UV. As shown 
in Table 1, the low ionization permitted lines, 
including O I A1304, C H A1335.3 and C H] 
AA2324.21,2328.84, Nil] AA2139.68,2143.45, AlII 
A1670.89 and A1 H] A2669.95, the various Si H 
multiplets, and Fe H 191 A1787 are all blueshifted 
by ^ 200 km s“^ with respect to z = 0.0608. 
Note that shifts of ^ 100 km s“^ are consis¬ 
tent with zero shift, given the accuracy of our 
wavelength zero point calibration using the in¬ 
terstellar medium lines. Higher ionization lines, 
including N HI] A1750, A1 HI AA1854.72,1862.79, 
Si HI] A1892.03 and Si IV AA1393.76,1402.77 are 
blueshifted by 300 — 500 km s while the 
highest ionization lines C IV AA1548.19,1550.77 
and N V AA1238.8,1242.8, are blueshifted by ~ 
900 km s“^. The highest blueshift, ~ 2000 km s“^, 
is displayed by He H A1640.4. 

This trend of increasing blueshift with increas¬ 
ing ionization level is observed in most quasars 
(e.g. Gaskell 1982; Wilkes 1984, 1986; Uomoto 
1984; Espey et al. 1989; Corbin 1990). How¬ 
ever, the typical amplitude of the high ioniza¬ 
tion line blueshift, such as C IV and N V, is 
only ~ 200 — 300 km s“^, and for He H it is 
~ 500 km s“^ (Tytler & Fan 1992; Laor et al. 
1995), i.e. about four times lower than found in 
I Zw 1. 

The low ionization line profiles are mostly 
consistent with the rather symmetric Ha pro¬ 
file. With increasing blueshift the lines get pro¬ 
gressively broader and develop a progressively 
stronger blue excess asymmetry (see Fig. 3). 


3.3. Reddening 

As shown in Fig.l, the optical-UV spectrum 
of I Zw 1 is significantly redder than observed in 
typical quasars. Could the intrinsic continuum of 
I Zw 1 be similar to the Francis et al. composite, 
and the observed red continuum just be due to 
reddening? 

Intrinsic reddening by Galactic like dust is 
ruled out as it would have produced a very 
pronounced broad absorption feature at 1800- 
2500 A, which is not observed. The grain size 
distribution in AGNs may, however, be different 
than in the Galactic ISM (e.g. Laor & Draine 
1993), and in particular it may not produce a 
significant 2200 A extinction bump. Thus, we 
cannot securely rule out reddening just based on 
the absence of the 2200 A feature. 

Reddening would also affect various optical/UV 
line ratios. In particular, assuming the correct 
optical/UV flux distribution is given by the Fran¬ 
cis et al. composite implies that Rq i ~ 1.2 (in¬ 
stead of 5.5), which is consistent with the value 
predicted by the Kwan & Krolik (1981) calcu¬ 
lation (which does not include the production 
of a O I] A1641 shown here to be insignificant). 
Dereddening also affects the Lya/H/? flux ratio. 
The intrinsic ratio would then be 32 instead of 
7.6. Such a high Lya/H/3 ratio is theoretically 
possible for plausible BLR parameters (e.g. Net- 
zer et al. 1995), but it is much higher than the 
ratio of 5-15 typically observed in quasars. 

If the intrinsic continuum of I Zw 1 is simi¬ 
lar to the Francis et al. composite, and the ob¬ 
served continuum is due to reddening by Galactic 
like dust (except the 2200 A bump) intrinsic to 
I Zw 1, then the implied reddening is of E(B- 
V) ~ 0.2. For a Galactic dust/gas ratio this im¬ 
plies Nu = 10^1 cm-A The ROSAT PSPG spec¬ 
trum of I Zw 1 implies an intrinsic neutral column 
density of no more than (1.5 ± 0.7) x 10^*^ cm”^ 
(see §4.2). Thus, the gas associated with the dust 
must either be highly ionized (U > 1) to prevent 
strong X-ray absorption below 0.4 keV, or the 
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dust/gas ratio needs to be at least 5 times larger 
than the Galactic value. The X-ray absorption 
by the grain component itself is negligible at the 
above columns (see Fig.6 in Laor & Draine 1993) 
and therefore cannot be used to constrain the 
dust column. 

A nearly featureless X-ray power-law may still 
be associated with a large dust column den¬ 
sity along the line of sight if most of the ob¬ 
served X-rays are electron scattered by an ex¬ 
tended medium. In this case only a small frac¬ 
tion of the X-rays would be scattered and thus al¬ 
though the X-ray power-law would be featureless, 
its normalization relative to the optical would 
be significantly reduced. However, Lawrence et 
ah (1997) find an optical to X-ray power-law of 
aox = —1.41, which is very close to the aver¬ 
age value of Uox = —1.48 for radio quiet quasars 
(Laor et ah 1997). Thus, it appears most likely 
that the primary X-ray source in I Zw 1 is ob¬ 
served directly, unless the direct optical flux is 
also completely obscured. 

Theoretical and observational uncertainties 
therefore do not allow us to clearly conclude 
whether I Zw 1 is reddened, or not. 

4. Intrinsic Absorption 

In this section we describe observational evi¬ 
dence for associated UV absorption and the im¬ 
plied constraints on the absorber column. We 
also discuss constraints based on the soft X-ray 
spectrum, and a possible relation of the absorber 
to a highly blueshifted emission component. 

4.1. UV absorption 
4-1.1. Observational evidence 

Figure 4 displays evidence for associated ab¬ 
sorption in Lyo, N V, and C IV. The clearest evi¬ 
dence for absorption is seen in the N V AA1238.8, 
1242.8 doublet. The N V absorption doublet is 
blueshifted by ~ 1870 km s“^ relative to the 
z = 0.0608 frame, it has a FWHM~ 300 km s“^, 
and an absorption EW of ~ 0.25 A (determined 


using the A1242.8 component). Since the spectral 
resolution of the FOS at N V is ~ 230 km s“^, 
the absorption system may not be truly resolved. 
The C IV AA1548.19,1550.77 profile is consistent 
with having a very similar absorption system, al¬ 
though the spectrum has a very low S/N at the 
blue wing of C IV (falls on the blue edge of the 
G190H). The Lya profile appears to be affected 
by an absorption system with the same velocity 
shift, and about the same EW and EWHM. 

The presence of high ionization lines, and lack 
of low ionization lines suggests the absorber is 
associated with I Zw 1, rather than with an in¬ 
tervening system unrelated to I Zw 1 (see e.g. 
Hamann 1997). Most of the constraints below 
are independent of the exact location of the ab¬ 
sorber. 

4-1.2. Implied constraints 

The column density associated with the ob¬ 
served absorption EW is a function of the ab¬ 
sorber optical depth. The optical depth can be 
deduced from the ratio of EW of the two doublet 
components, even when the absorption line pro- 
hles are not resolved, assuming the absorber com¬ 
pletely covers the emission source. In the opti¬ 
cally thin limit this ratio equals the ratio of oscil¬ 
lator strengths, i.e EW(1238.8)/EW(1242.8)= 2, 
and when both absorption lines are saturated 
this ratio approaches unity. The observed ra¬ 
tio is ~ 1.4, but it is rather uncertain since the 
1242.8 A component absorbs the line peak where 
the underlying emission is not clearly defined, 
and the 1238.8 A component is partly blended 
with Galactic Si II A1304.37 absorption. The ab¬ 
sorption optical depth therefore remains uncer¬ 
tain. 

In the optically thin limit the absorber column 
density is related to the EW through 

iVabs = 1.13 X 102°A-V^EW cm-A 

where A, the absorption wavelength, and its EW 
are in units of A, and / is the oscillator strength. 
Using /(1215.67) = 0.416,/(1550.77) = 0.095, 
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/(1242.8) = 0.078 and an absorption EW of 
0.25 A for all lines, we get A^h i = 4.6x 10 ^^ cm“^, 
IV = 1-2 X 10^^ cm“^, and A^n v = 2.3 x 
10^^ cm“^. These values are lower limits on the 
absorber columns. 

The upper limits on the columns are obtained 
in the optically thick limit, i.e. when only ther¬ 
mal broadening is present. If the absorption line 
profile is not resolved, and only the absorption 
EW is given, then the deduced column density 
depends on whether the line is optically thick in 
the Doppler core only, or whether it is also thick 
in the Lorentzian wings (i.e. “damped absorp¬ 
tion”). The first case applies when 

EW/AXd ^ 6 , 

where 

AXd = 4.285 X 

is the thermal line width in A, T = 10 ^T 4 K, and 
A is the atomic weight. The absorber column 
density is then given by 

A^abs = 2 X 

When absorption is dominated by the Lorentzian 
wings, i.e. when 

EW/AXd 6 , 

The absorber column density is 

A^abs - 3.35 X 103®A“^/“^r-^EW^ 

where T is the radiative decay rate (ri 2 i 5.67 = 
6.265x108 S-l,ri550.77 = 2.64x108 S-l,ri242.80 = 
3.387 X 108 Assuming T 4 = 2, the typi¬ 

cal temperature of photoionized gas where C IV 
or N V dominate, we get the maximum possi¬ 
ble columns of Vh i = 4.3 x 10^^ cm“^ (where 
EW/AA^) < 6 ) Nq IV = 1.4 X 10 ^® cm“^, and 
A^n V = 3.3 X 10^8 cm“^ (both cases EW/AAd > 
6 ). The H I column density upper limit is rather 
uncertain since a Doppler core absorption gives 
Aabs oc and thus a likely ±30% 


error in the H I EW corresponds to an upper 
limit in the range of 0.1 — 3.1 x 10^® cm“^. The 
C IV and N V upper limits are uncertain by only 
a factor of 2 , since Vabs oc EW^, when the 
Lorentzian wings dominate the absorption. 

To infer the total H column density from the 
H I column density one needs to know the H ion¬ 
ization state. The H II/H I fraction in the ab¬ 
sorbing gas is related to the ionization parameter 
through 

A^h ii/Ah I — 2 X 10^1/ 

(e.g. Netzer 1990). Thus, the upper limit on the 
H column density is 

Vh — A^h II = (0.2 — 6 ) X 10^°[/ cm“^. 

The presence of N V and C IV ions suggests that 
the absorber has U ~ 0.01 — 1, thus the H col¬ 
umn density upper limit is probably in the range 
10^^ ^ Ah ^ 6 X 10^°. If the absorber metal 
abundance is solar then the above Ah implies 
metal columns upper limits of 3.6 x 10 ^^ ^ jVc ^ 
2 X 10^^ and lO^^ ^ ^ 6.7 x 10^% These 

values overlap with the directly determined con¬ 
straints on the C IV and N V columns and sug¬ 
gests that significant fractions of C and N are 
indeed in the form of C IV and N V, consistent 
with the assumption of 1/ ~ 0.01 — 1. Detection 
of O VI AA1031.93,1037.63 absorption would al¬ 
low a much tighter constraint on the absorber’s 
U. 

4.2. Soft X-ray absorption 

In addition to the UV absorption lines there 
may also be some soft X-ray absorption in I Zw 1. 
Boiler, Brandt, & Eink (1996) and Lawrence et 
al. (1997) fit a single power-law to the ROSAT 
PSPC spectrum of I Zw 1 and find a best fit 
column density of Ah i = (6.5 — 6.7 ± 0.7) x 
10^® cm“^ (90% error) versus a Galactic value of 
(5.0 ± 0.1) X 10^° cm“^ (Elvis et al. 1989), sug¬ 
gesting that there may be an intrinsic absorber 
in I Zw 1 with Ah ~ (1.5 —1.7±0.7) x 10^° cm“^. 
Lawrence et al. (1997) find that a fit including 
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only Galactic absorption can be excluded with a 
confidence level greater than 99%. 

The observed total column density in I Zw 1 
is optically thick below 0.38 keV, where He dom¬ 
inates the opacity. Thus, the only direct con¬ 
straint from ROSAT is for iVue i ~ (1.5 ± 0.7) x 
10 ^® cm“^, assuming the cosmic = 0.1 

abundance. A ~ 20% higher column density ap¬ 
plies if most of He is He H, as deduced using the 
Verner & Yakovlev (1995) He I and He H absorp¬ 
tion cross sections. 

Associations of UV and X-ray absorbers with 
the same gas was claimed in 3C 351 (Mathur 
et al. 1994), and NGC 5548 (Mathur, Elvis & 

Wilkes 1995), but disputed in NGG 3516 (Kriss 
et al. 1996a, 1996b). Gan the UV and the X- 
ray absorption in I Zw 1 originate in the same 
absorber? In a cloud photoionized by a typical 
AGN continuum there is a surface layer of thick¬ 
ness Vh ~ 10^^ [/ cm“^ where H is mostly ion¬ 
ized. At the surface of this “H H” layer there 
is a thinner “He HI” layer. For a typical AGN 
continuum the column density of the “He HI” 
layer is VH(He HI) ~ 2 x 10^^17 cm“^. Thus, 
irrespective of the value of U, the “He HI” layer 
is about three times thicker than the upper limit 
of Nil < 6 X lO^^t/ cm“^ deduced from the UV 
absorber in I Zw 1. The He in the UV ab¬ 
sorber would be mostly He HI, and the soft X- 
ray absorption would originate in the residual 
He H within the UV absorber. The He Hi/He H 
fraction within the “He HI” region is given by 
A^He iii/A^He II — 10^1/ (Netzer 1990). The up¬ 
per limit on the He H column density in the UV 
absorber is then 

NRe II = 1VhX^x^%^ < 6x102°[/x0.1x(103U)“\ 

^ Nn Nne ~ y J ^ 

i.e. Nile II = 6x 10^® cm“^ which is ~ 300 smaller 
than the Nue ii required to produce the X-ray 
absorption. It thus appear that the upper limit 
on the column of the layer which produces the 
UV absorption is much smaller than required to 
produce the possible X-ray absorption. 


This line of argument suggests that the UV 
and X-ray absorbers are distinct. However, given 
the uncertainty in the upper limit on the Lya 
absorption EW, which sets Nn, and the uncer¬ 
tainty in the far UV spectral shape, which sets 
A^He iii(U), we cannot conclusively rule out the 
possibility that an “He H” region is present be¬ 
hind the “He HI” region. Such an “He H” region 
can produce a high He H column, which would be 
required if the X-ray absorption is real, together 
with a relatively small total H I column, implied 
by the UV absorption lines. 

4.3. UV and optical emission 

The UV absorber blueshift of 1870 km s ^ is 
remarkably close to the He H A1640 emission line 
blueshift of ~ 1990 km s“^, suggesting that the 
strongly blueshifted He H emission line peak may 
originate in the same outflowing gas which pro¬ 
duces the UV absorption lines. In order to emit 
the observed He H A1640 flux the absorber must 
absorb most photons just above the He H bound- 
free edge at 4 Rydberg. As described above, 
when Nn ~ 6 x 10^^U cm“^ the UV absorber 
is within a factor of three of the column density 
required to form an “He H” region, i.e. it absorbs 
~ 1/3 of He H ionizing photons, and thus if the 
column density of the UV absorber is close to 
the upper limit, and its covering factor is close 
to unity, it may contribute significantly to the 
observed He H A1640 emission. 

The intrinsic X-ray absorber becomes opti¬ 
cally thick below 0.2 keV. It is therefore highly 
optically thick at 4 Rydberg and can produce the 
observed He H A1640 emission if its covering fac¬ 
tor is reasonably large (> 10%). 

It is interesting to note that Oke & Lauer 
(1979) found evidence for [O HI] A5007 and 
[Ne HI] A3868 emission components at z = 0.0548 
which corresponds to a 1800 km s ^ blueshift 
with respect to the z = 0.0608 used here. The 
same feature was also noted by van Groningen 
(1993) in a high spectral resolution profile of 
[O HI] A5007. If these emission components also 
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originate in the UV absorber then their flux can 
be used as an additional constraint on the den¬ 
sity, column density, and ionization parameter of 
the UV absorber. 

4.4. Future observations 

Higher spectral resolution observations of the 
Lya N V and C IV absorption lines would al¬ 
low a much more accurate determination of the 
absorption profile, and thus of the absorber’s col¬ 
umn. If the H I column density is indeed close to 
the upper limit, as required to “unite” the UV 
absorption, X-ray absorption, and He H emis¬ 
sion, then the Lya absorption profile, for exam¬ 
ple, would appear as a ~ 60 km s“^ wide, steep 
sided, and flat bottomed trough going down to 
zero flux. The metal absorption lines should also 
be resolved as deep narrow optically thick ab¬ 
sorption lines with damping wings. If the ab¬ 
sorption lines remain shallow and broad then the 
lower limits to the absorbing columns would ap¬ 
ply. The metal lines absorption doublet ratios 
can be used to establish the covering factor of 
the UV absorber. 

Higher S/N X-ray spectra below 0.5 keV are 
required to establish whether the soft X-ray ab¬ 
sorption hinted at by the i^O^HTPSPC is indeed 
real. 

5. Discussion 

The strong optical Fe H emission of I Zw 1, its 
weak [O HI] emission, strong IR emission, and 
“red” UV continuum are all typical properties 
of low ionization broad absorption line quasars 
(BALQSOs; Weymann et al. 1991, Boroson & 
Meyers 1992; Sprayberry & Foltz 1992). The 
presence of weak UV absorption at a blueshift 
of ~ 2000 km s“^ in I Zw 1 suggests it may be a 
“failed” low ionization BALQSO, i.e. a BALQSO 
where our line of sight just grazes the outflow¬ 
ing high ionization wind, and misses the low ion¬ 
ization outflow, as suggested by Turnshek et al. 
(1994) in the case of PG 0043-1-039. 


Is the red continuum, excess blue flux in the 
high ionization lines, associated absorption, and 
dense BLR, unique to I Zw 1, or are they typical 
UV properties of narrow line quasars? No com¬ 
plete UV study of narrow line AGNs is available 
to answer these questions, but some hints may 
be obtained from existing observations (see also 
discussion in Lawrence et al. 1997). 

Baldwin et al. (1996) studied the emission 
line properties of a heterogeneous sample of seven 
z ~ 2 quasars with a range of emission line prop¬ 
erties. Some of their objects have strong Fe H, 
Fe HI, and Al HI emission as observed in I Zw 1. 
Baldwin et al. analyzed in detail the UV spec¬ 
trum of Q0207 — 398, where a high S/N ratio 
was available, and found excess emission in the 
blue wing of the high ionization UV lines, which 
they interpreted as a high ionization outflowing 
component in the BLR, although there was no 
direct evidence (through absorption) for such a 
component, as probably observed in I Zw 1. 

It is also interesting to note that BG92 found 
that strong optical Fe H emitting quasars, which 
tend to have narrow H/3, also tend to have blue 
excess flux in H/3, and Boroson & Meyers found 
that low ionization BALQSOs, which are gener¬ 
ally I Zw 1 like, have blue excess flux in Ha. It 
is not known, however, whether this property ex¬ 
tends to the UV lines as well. 

If the UV absorption system in I Zw I is in¬ 
deed producing the blueshifted components of 
the high ionization lines, then the absence of 
a corresponding redshifted component indicates 
that the far side of this outflow has to be ob¬ 
scured. This may either be due to an extended 
highly optically thick gas, such as an accretion 
disk, or it may result from absorption within 
each cloud, if the clouds column density is large 
enough (see Ferland et al. 1992 for discussion). 

The doublet ratios observed here indicate that 
the low ionization lines are thermalized and will 
not be emitted isotropically. However, the large 
velocity gradients in the wind may strongly re¬ 
duce the optical depth for the high ionization 


16 



lines in the wind. The observed blueshifted emis¬ 
sion in C IV and N V is far too broad to deter¬ 
mine the observed doublet ratio and cannot be 
used to determine if the emission from that gas 
is optically thick or thin. 

I Zw 1 has particularly weak C III] A 1909 
emission (§3.1.2) but normal Si III] A 1892 emis¬ 
sion, which we interpret here as evidence for a 
relatively dense ( 10^^ cm BLR, at least for 

the region which produces the the low ionization 
lines. Baldwin et al. (1988) noted the weakness 
of C III] A 1909 in four other quasars with nar¬ 
row UV lines, and this trend is also apparent in 
Baldwin et al. (1996) quasars. The most ex¬ 
treme case is of H0335 — 336 where essentially 
no C III] A 1909 emission was observed (Hartig 
& Baldwin 1986). This quasar has very narrow 
lines, and is also a low ionization BALQSO. 

Independent evidence for a dense BLR with 
Ue ^ 10^^ cm“^ is provided by the significant EW 
and thermalized doublet ratio of the Al III A1857 
doublet (§3.1.6), and by the significant EW of 
C IIP A1176 (§3.1.2). Another indication for a 
dense BLR in I Zw 1 is provided by its optical 
spectrum. The Na I AA5889.95, 5895.92 emission 
line is rather strong (Oke & Lauer 1979; Phillips 
1976), which Thompson (1991) and Korista et al. 
(1997, their Eigure 3g) find can only be produced 
for Ue ^ 10^^ cm“^, and a large column density 
(required to shield neutral Na I from ionizing ra¬ 
diation at E > 5.14 eV). 

I Zw 1 has a “red” continuum (§3.3), which 
is rather rare among quasars. Two of the four 
quasars with narrow lines described by Baldwin 
et al. 1988 also have a a — 1 UV spectral slope, 
but the other two narrow line quasars have a ~ 0 
and —0.5, which are typical values. 

It thus appears that some of the properties 
of I Zw 1, in particular the relatively weak 
C III] A 1909, and the blueshifted excess emission 
in the high ionization lines, may be common in 
narrow line quasars. A more systematic study of 
the UV emission of narrow line AGNs is required 
to establish their typical emission line properties. 


their relation to low ionization BALQSOs, and to 
eventually understand the underlying physics. 

6. Summary 

A high S/N UV spectrum of I Zw 1 was ob¬ 
tained with the HST EOS with a spectral res¬ 
olution of A/AA ~ 1300 over the 1150-3250 A 
spectral region. The following main results are 
obtained from our analysis of the UV spectrum 
of I Zw 1. 

1. The Mg II A2798 doublet is partially re¬ 
solved. The measured doublet ratio is 
1 .2/1, consistent with theoretical predic¬ 
tions that the BLR is highly optically thick 
to absorption in this line. 

2. The Al III A1857 doublet is clearly resolved 
with an “optically thick” doublet ratio of 
1.25/1. The line optical depth provides an 
upper limit to the clouds distance from the 
ionizing source which is consistent with the 
“standard” BLR radius. The strength of 
this line requires a BLR component with 
rie ^ 4 X 10^^ — 4 X 10^^ cm“^. 

3. A weak UV absorption system is detected 
in N V, and possibly C IV and Lya as well, 
indicating an outflow with a line of sight 
velocity of ~ 1870 km s“^ and velocity dis¬ 
persion 300 km s“^. 

4. Lines from ions of increasing ionization 
level show increasing excess flux on the blue 
wing, and an increasing shift in the veloc¬ 
ity of the emission line peak, reaching a 
maximum blueshift of 2000 km s ^ for 
He II A1640. This may indicate an outflow¬ 
ing component in the BLR, visible only in 
the approaching direction, where the ion¬ 
ization level increases with velocity. The 
highest velocity part of this outflow may 
produce the observed UV absorption sys¬ 
tem. 
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5. The small C III] A1909 EW, and the small 

C III] AI909/Lya and C III] A1909/Si III] A1892 
flux ratios indicate that the bulk of the gas 
in the BLR has a density of 10^^ cm“^, 
about an order of magnitude larger than 
typically seen in AGNs. 

6 . Very weak [C III] A1907 and [Si III] A1883 
emission may be present, indicating there 
may be a NLR component with ne ~ 5 x 
10^ cm“^. 

7. The C III* A1176 line is clearly detected. 

Its strength indicates a BLR component 
with Ue ^ 10^^ cm“^. 

8 . Prominent Fe II UV 191 emission is seen, 
together with weaker emission lines at 1294 A 
and 1871 A. These three features have been 
proposed in the literature as evidence for 
significant Lya pumping of the 8-10 eV lev¬ 
els of Fe 11. 

9. Significant Fe III emission is present. Fe III 
UV 34 and UV 48 multiplets are clearly re¬ 
solved, and Fe III UV 1, UV 47, UV 50, 
and UV 68 may also be present. The im¬ 
plications of significant Fe III emission for 
the conditions in the BLR needs to be ex¬ 
plored. 

10. The suggestion by Grandi (1983) that the 
O I A1304 emission line is partly converted 
to O I] A1641 is ruled out in the case of 
I Zw 1, indicating that the bulk of O I 
Ly/I Bowen flouresence occurs in the BLR 
clouds at a depth where TLya < 10®. The 
strong suppression of O I A1304 in I Zw 1 is 
either due to reddening, or to Balmer con¬ 
tinuum absorption. 

The HST spectrum of I Zw 1 demonstrates 
the richness and complexity of quasars emission 
line spectra. It is not clear whether some of 
the emission features are unique to I Zw 1, or 
whether they are typical of most AGNs. Many of 
the features detected here would be blended and 


unidentifiable in typical quasar spectra, creating 
a deceptively simple emission spectrum. A close 
inspection of the Francis et al. composite (Fig.l) 
suggests that most of the features detected here 
are indeed present in other AGNs. The rich UV 
emission spectrum of I Zw 1 should serve as a 
useful benchmark for photoionization codes, in 
particular for codes attempting to calculate the 
complex emission spectrum of Fe 11. Many weak 
features remain unidentified, and it is likely that 
at least some of these are weak Fe II features. 
A reliable identification of these weak features 
may follow from detailed photoionization models. 
The high quality optical-UV spectrum of I Zw 1 
obtained in this study, can be used for construct¬ 
ing an Fe II UV template, and is available upon 
request from the authors. 
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Table 1: Emission Line Parameters^ 


■^obs 

A 

^ID 

A 

Ion 

Flux'’ 

EW 

A 

FWHM 
km s“^ 

Ay 

km s“^ 

Comments 

1110-1130 

1122-I-? 

Fe III UV 14-? 

31.1 

12.6 




1174.9 

1175.7 

C III* 

5.2 

1.8 

1210 

-200 


1206.2 

1206.5 

Si III 




-75 

affected by absorption on both sides 

1214.2 

1215.7 

H I 

382 

128 

1460 

-370 

N V subtraction is only approximate 

1236.7, 1239.1 

1238.8, 1242.8 

N V 

52 

19 

2350 

-890 

velocity shift based on 1242.SA component 

1249.5 

1246.7, 1248.4, 1251.1 

Si II* 





tentative identification 

1260.5, 1264.5 

1260.4, 1264.7, 1265.0 

Si II 

10.5 

3.1 


0 

strongly blended components 

1302.6, 1305.5 

1302.2, 1304.9, 1306.0 

0 I 

20.6 

8.3 


0 

flux and EW for O I -1- Si II blend 

1309.2 

1304.4, 1309.3 

Si II 

20.6 

8.3 


0 

as above 

1334.3 

1335.3 

C II 

8.8 

3.6 

1250 

-220 

unidentified 

1344.3 

7 

7 

2.2 

0.9 

1050 



1392.1, 1400.6 

1393.76, 1402.77 

Si IV 

47.7 

17.4 


-360 

flux and EW for Si IV + O IV] blend 


1401 

0 IV 

47.7 

17.4 



as above 

1543.4, 1546.0 

1548.19, 1550.77 

C IV 

107 

28.9 

3930 

-930 


1629.5 

1640.4 

He II 

9 

2.5 

1700 

-1990 

blended features at 1620A and 1640A ? 


1660.81, 1666.15 

O III] 

8.3 

2.5 



flux and EW for O III] -|- A1 II blend 

1670.5 

1670.89 

A1II 

8.3 

2.5 


-70 

flux and EW for O III] + A1 II blend 

1746.9 

1749.7 

N III] 

12.5 

3.5 

1800 

-480 

broad multiplet (1746.82A-1753.99A) 

1785.8 

1785.3, 1786.7, 1788.0 

Fe II UV 191 

4.8 

1.4 

910 

-150 

approximate shift 

1807.9, 1816.2 

1808.01, 1816.93, 1817.45 

Si II 

5.0 

1.5 


0 


1852.7, 1860.8 

1854.72, 1862.79 

A1 III 

16 

4.4 


-320 

blended features at 1840A and 1870A 

1881.1 

1882.7 

[Si III]? 




-250 

marginal detection 

1890.9 

1892.03 

Si III] 

32.5 

9.1 


-250 

blended with C III] 

1894.0 

1895.46 

Fe III UV 34 




-230 

marginal detection 

1906.4 

1906.68 

[C III]? 




-40 

as above 

1907.8 

1908.73 

C III] 

17.7 

5.1 


-150 

flux and EW for the C III]-FFe III UV 34 1914 blend 

1912.6 

1914.06 

Fe III UV 34 

17.7 

5.1 


-230 

as above 

1924.4 

1926.30 

Fe III UV 34 

7.6 

2.2 


-300 


1951.7 

1951.0, 1952.6, 1953.3 

Fe III UV 68 

5.7 

1.8 


-140 

shift relative to central component 

2060.3, 2067.2 

2062.20, 2068.90 

Fe III UV 48 

6.2 

2.0 


-260 


2078.9 

2079.65 

Fe III UV 48 

1.4 

0.5 

640 

-110 


2141.4 

2139.68, 2143.45 

N II] 

6.5 
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Fig. 1.— The optical-UV emission of I Zw 1, corrected for Galactic reddening (for E(B—V)=0.1), versus 
the Francis et al. composite (set to match I Zw 1 at logi^ = 14.7). Note that there is practically no 
wavelength range without significant emission features. Many of the very weak features in the Francis et 
al. composite are clearly resolved in the I Zw 1 spectrum. The optical spectrum of I Zw 1 was available 
through the courtesy of B. Wills. 
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Fig. 2.— a. Tentative identifications of most of the emission features in the HST FOS spectrum of 
I Zw 1. The lines indicate the expected wavelength for z = 0.0608. Only a small number of the Fe II 
multiplets which may be present were marked. Line designations below the spectrum refer to Galactic 
absorption features. 
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Fig. 2.— d. The optical spectrum. Note that the absolute flux calibration is not secure. The overlapping 
HST spectrum is plotted as a dotted line at A < 3100 A. This spectrum was shifted down by 0.37 in log 
flux to match the optical spectrum. 
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Fig. 3.— Emission line profiles for selected lines observed in the spectrum of I Zw 1. The histogram 
represents the observed profile and the solid line is the Ha profile fit. Multiplets are ht with a sum of 
template components. The wavelength used to set the velocity scale and the multiplet component flux 
ratios are indicated in each panel. Note the overall bad fit and large velocity shifts of the high ionization 
lines (upper panels), and the much better fit and very small shifts of the low ionization lines (lower 
panels). 
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Fig. 4.— Absorption line profiles of the candidate intrinsic UV absorption line system. The rest 
wavelength used to set the velocity scale is indicated in each panel. The strongest evidence for absorption 
is seen in the N V A1242.80 profile. The other lines are consistent with having absorption of similar EW 
and FWHM at the same velocity shift of ~ 1870 km s“^. 
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